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Abstract 

The quantum interference effects induced by the Wess-Zumino term, or Berry 
phase are studied theoretically in resonant quantum coherence of the mag- 
netization vector between degenerate states in nanometer-scale single-domain 
ferromagnets in the absence of an external magnetic field. We consider the 
magnetocrystalline anisotropy with trigonal, tetragonal and hexagonal crys- 
tal symmetry, respectively. By applying the periodic instanton method in the 
spin-coherent-state path integral, we evaluate the low-lying tunnel splittings 
between degenerate excited states of neighboring wells. And the low-lying 
energy level spectrum of m-th excited state are obtained with the help of 
the Bloch theorem in one-dimensional periodic potential. The energy level 
spectrum and the thermodynamic properties of magnetic tunneling states are 
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found to depend significantly on the total spins of ferromagnets at sufficiently 
low temperatures. Possible relevance to experiments is also discussed. 
PACS number(s): 75.45.+j, 75.10.Jm, 03.65.Bz 
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I. INTRODUCTION 



In recent years there has been great experimental and theoretical effort to observe and 
interpret macroscopic quantum tunneling (MQT) and coherence (MQC) in nanometer-scale 
single-domain magnets. One notable subject is that the topological Berry or Wess-Zumino 
phase&i can lead to remarkable spin-parity effects. Loss et a/.,@ and von Delft and Henlyi 
showed that the tunnel splitting is suppressed to zero for half-integer total spins in biaxial 
ferromagnetic (FM) particles due to the destructive phase interfererence between topologi- 
cally different tunneling paths. However, the phase interference is constructive for integer 
spins, and hence the splitting is nonzerofii While spin-parity effects are sometimes be re- 
lated to Kramers degeneracyjl! they typically go beyond the Kramers theorem in a rather 
unexpected way.§i Barnes et al. proposed the auxiliary particle method to study the model 
for a single large spin subject to the external and anisotropy fields, and discussed the spin- 
parity effects.! Similar effect was found in antiferromagnetic (AFM) particles, where only 
the integer excess spins can tunnel but not the half-integer ones.0lll Recently, topological 
phase interference effects were investigated extensively in FM and AFM particles in a mag- 
netic field,i'i0lii'0 and in the systems with different symmetriesJilHll Spin tunneling and 
quantum oscillation at excited states were studied for biaxial FM particles at zero mag- 
netic field,0 and at a field along the hard axisS One recent experiment^ was performed to 
measure the tunnel splittings in molecules Fes, and a clear oscillation of the splitting as a 
function of the field along the hard axis was observed, which is a direct evidence of the role 
of the topological spin phase (Berry phase) in the spin dynamics of these molecules. 

It is noted that the previous results of topological phase interference effects were ob- 
tained for the tunnel splittings of the ground state in FM particles with different crys- 
tal symmetries, El or for the excited states in FM particles with simple biaxial crystal 
symmetry.0 The purpose of this paper is to study the spin-parity effects at excited states for 
FM particles with a more complex (than biaxial) structure, such as trigonal, tetragonal, and 
hexagonal symmetry around z, which have three, four, and six degenerate easy directions 
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in the basal plane. Integrating out the momentum in the path integral, the spin tunneling 
problem is mapped onto a particle moving problem in one-dimensional periodic potential 
V (4>). By applying the periodic instanton method, we obtain the low-lying tunnel splittings 
between m-th degenerate excited states of neighboring wells. The periodic potential V ((f)) 
can be regarded as a one-dimensional superlattice. The general translation symmetry results 
in the energy band structure, and the low-lying energy level spectrum of excited states is 
obtained by using the Bloch theorem and the tight-binding approximation. Our results show 
that the excited-state tunnel splittings depend significantly on the parity of the total spins. 
And the structure of energy level spectrum for the trigonal, tetragonal and hexagonal crystal 
symmetry is found to be much more complex than that for the biaxial crystal symmetry. 
Another important conclusion is that the spin-parity effects can be reflected in thermody- 
namic quantities of the low-lying tunneling levels. Thermodynamic property (such as the 
specific heat) of the magnetic tunneling states is evaluated, and is found to be strongly par- 
ity dependent on the total spins, which may provide an experimental test for the topological 
phase interference effects. And the spin-parity effect is lost at high temperatures. 

The remaining part of this paper is organized as follows. In Sec. II, we review briefly some 
basic ideas of MQT and MQC in FM particles, and discuss the fundamentals concerning 
the computation of excited-level splittings in the double-well-like potential. In Sees. Ill and 
IV, we study the spin tunneling between degenerate excited states in FM particles with the 
trigonal, tetragonal and hexagonal symmetry. The conclusions are presented in Sec. V. 



II. SPIN TUNNELING IN FM PARTICLES 

For a spin tunneling problem, the tunnel splitting for MQC or the decay rate for MQT is 
determined by the imaginary-time transition amplitude from an initial state \i) at r = —T/2 
to a final state |/) at r = T/2 in the spin-coherent-state representation aJE>0 

U fi = (f\ e~ HT \i) = J Dfiexp {-S E ) , (1) 
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where DQ = sm9d9d(f>. The paths appearing in Eq. (1) are fixed at the end points r = 
±T/2. For a system with equivalent double wells, we let \E) + and \E)_ be eigenstates of 
the same energy E in the right- and left-hand wells, respectively. The small contribution 
due to quantum tunneling leads to the effect of level splitting AE, which removes the 
asymptotic degeneracy. The corresponding eigenstates of the Hamiltonian separate into 
odd and even states \E) and \E) e which are superpositions oi\E) + , \E)_ such that \E) — 
^5 (\E) + - \E)_), and \E) e = ^ (\E) + + \E)_) with eigenvalues E ± AE, respectively. 
In the limit that T — > oo, one expects that the amplitude for the transition from state 
\E)_ in the left-hand well to the state \E) + in the right-hand well in the time interval T as 
+ (E\ e~ HT \E)_ — > exp (-ET) sinh (AET). Therefore the tunnel splitting AE is obtained if 
the transition amplitude can be calculated. The decay rate T from the metastable state for 
MQT can be evaluated by a similar procedure. The Euclidean action Se in Eq. (1) is!i3§! 



S E (6,<P) = j I ,It 



.Mo (d<f>\ .M d<p\ 
i — I — I - % — — cos + E 
7 Vary 7 Vary 



(2) 



where M = M — h'jS/V, V is the volume of the particle, 7 is the gyromagnetic ratio, 
and S is the total spins. It is noted that the first two terms in Eq. (2) define the Berry 
or Wess-Zumino term which has a simple topological interpretation. For a closed path, this 
term equals —iS times the area swept out on the unit sphere between the path and the 
north pole. The first term in Eq. (2) is a total imaginary-time derivative, which has no 
effect on the classical equations of motion, but it is crucial for the spin-parity effects. i@0 

In the semiclassical limit, the instanton's contribution to T or AE (not including the 
topological Wess-Zumino phase) is given by0 

T(or AE)=Au p (^j 2 e- s °\ (3) 

where u p is the oscillation frequency in the well, S c i is the classical action, and the prefactor 
A originates from the quantum fluctuations about the classical path. It is noted that Eq. 
(3) is based on tunneling at the ground state, and the temperature dependence of the tun- 
neling frequency (i.e., tunneling at excited states) is not taken into account. The instanton 



technique is suitable only for the evaluation of the tunneling rate at the vacuum level, since 
the usual (vacuum) instantons satisfy the vacuum boundary conditions. Recently, Liang et 
aZ.il developed new types of pseudoparticle configurations which satisfy periodic boundary 
condition (i.e., periodic instantons or nonvacuum instantons). For a particle moving in a 
double- well-like potential U (x) , the WKB method gives the tunnel splitting of excited states 
at an energy E > aS§H! 

AE = ^1 exp [-S (E)] , (4) 

7T 

with the imaginary-time action is 

rx 2 (E) i 

S(E) = 2v / 2^/ dxJU(x)-E, (5) 

Jxi(E) v 

where {E) are the turning points for the particle oscillating in the inverted potential 
— U (x). to (E) = 2n/t (E) is the energy- dependent frequency, and t (E) is the period of the 
real-time oscillation in the potential well, 

/-— f x *( E ) dx , s 

t (E) = / 6 

Jx 3 (E) JE-U (x) 



where X3 4 (E) are the classical turning points for the particle oscillating inside U (x). The 
functional-integral and the WKB method showed that for the potentials parabolic near the 
bottom the result Eq. (4) should be multiplied by J%^^^\~~ S'S This factor is very close 
to 1 for all n: 1.075 for n = 0, 1.028 for n = 1, 1.017 for n = 2, etc. Stirling's formula for 
n! shows that this factor trends to 1 as n — ► 00. Therefore, this correction factor, however, 
does not change much in front of the exponentially small action term in Eq. (4). Recently, 
the crossover from quantum to classical behavior and the associated phase transition have 
been investigated extensively in nanospin systems&^ and other systems.0 

III. MQC FOR TRIGONAL CRYSTAL SYMMETRY 



In this section, we consider a spin system with trigonal crystal symmetry, i.e., which has 
three consecutive energy minima in a period. Now the total energy is 
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E (9, 0) = K x cos 2 9 - K 2 sin 3 9 cos (30) + E , (7) 

where K 1 and K 2 are the magnetic anisotropic constants satisfying K\ ^> K 2 > 0, and 
E is a constant which makes E (9, 0) zero at the initial state. As K x 3> K 2 > 0, the 
magnetization vector is forced to lie in the 9 = n/2 plane, so the fluctuations of 9 about n/2 
are small. Introducing 9 = n /2 + a (\a\ <^ 1), Eq. (7) reduces to 

E (a, 0) « K x a 2 + 2X 2 sin 2 (30/2) . (8) 

The ground state corresponds to the magnetization vector pointing in one of the three 
degenerate easy directions: 9 = ir/2, and = 0, 27r/3, 47r/3, other energy minima repeat 
the three states with period 2n. Performing the Gaussian integration over a, we can map 
the spin system onto a particle moving problem in one-dimensional potential well. Now the 
transition amplitude becomes 

U f i = exp [-iS (0/ - 4>i)\ j d0exp (S E [0]) , 
= exp [—iS((f>f — (f>i)] J <i0exp|— J dr 

with m = hS 2 /2K 1 V, and V(<j>)=2 (K 2 V/k) sin 2 (30/2). It is noted that the total deriva- 
tive in Eq. (2), when integrated, gives an additional phase factor to the transition amplitude 
Eq. (9) which depends on the initial and final values of 0. For the trigonal symmetry, this 
phase factor in Eq. (9) is exp (—i2nS/?>). The potential V (0) is periodic with period 27r/3, 
and there are three minima in the entire region 27r. We may look at V (0) as a superlattice 
with lattice constant 27r/3 and total length 2n, and we can derive the energy spectrum by ap- 
plying the Bloch theorem and the tight-binding approximation. The translational symmetry 
is ensured by the possibility of successive 2tt extensions. 

The periodic instanton configuration P which minimizes the Euclidean action in Eq. (9) 
satisfies the equation of motion 

\m(^j -V = (10) 
6 



-m 



dr 



(9) 



where E > is a constant of integration, which can be viewed as the classical energy of the 
pseudoparticle configuration. Then we obtain the kink-solution as 



si ir [ -(j) p ) = 1 — A; 2 sn 2 (uit, k) , 



(11) 



where sn^r, k) is the Jacobian elliptic sine function of modulus k, 



p < - 1 

o 1 

nf 



(12) 



with u! = 3y/2(V/hS)y/K^, and m = ^2K 2 V/hE > 1. In the low energy limit, i.e., 
E — > 0, k — > 1, sn(w, 1) — > tanh-u, we have 

' 3 ' ^ 1 (13) 



S1 " W'J - cosh 2 ( Wl r) ' 
which is exactly the vacuum instanton solution derived in Ref. 16. 

The Euclidean action of the periodic instanton configuration Eq. (11) over the domain 
(— 13,13) is found to be 



S p — 



dr 



1 d<p 

2 I dr 



p 



= W + 2E(3, 



(14) 



with 



W 



2 5/2 q [K 2 
3 



E{k) - (l-k 2 )K(k) 



(15) 



where K (k) and E (k) are the complete elliptic integral of the first and second kind, respec- 
tively. Now we discuss the low energy limit where E is much less than the barrier height. 
In this case, k' 2 = 1 — k 2 = hE/2K 2 V < 1, so we can perform the expansions of K (k) and 
E(k) in Eq. (15) to include terms like k! 2 and k! 2 In (4/A;'), 



E(k) =1 + - 



h 4 



r 

2. 



k' 2 + 



k' 2 + 



(16) 



With the help of small oscillator approximation for energy near the bottom of the potential 
well, E = e l ™ = (m + 1/2) Eq. (15) is expanded as 



25/2 fj£ , y\ / \ 

W = —J-S-(m + -) + (m + -)ln 



IK! If 1 
m + 



2 9 / 2 y KoS\ 2 



(17) 



Then the general formula Eq. (4) gives the low-lying energy shift of m-th excited states for 
FM particles with trigonal crystal symmetry at zero magnetic field as 

h*e% = -L ; (^V) A*Ki) 5 (-*) exp (~X 1/2 s) , (18) 

where A = K 2 /Ki, and 5 is the total spin. 

It is noted that hAe 1 ™ is only the level shift induced by tunneling between degenerate 
excited states through a single barrier. The periodic potential V (0) can be regarded as a 
one-dimensional superlattice. The general translation symmetry results in the energy band 
structure, and the energy spectrum could be determined by the Bloch theorem. It is easy 
to show that if e l ™ are the degenerate eigenvalues of the system with infinitely high barrier, 
the energy level spectrum is given by the following formula with the help of tight-binding 
approximation, 

E™ = e>™ - 2Ae« cos [(S + 2vr/3] . (19) 

The Bloch wave vector £ can be assumed to take either of the three values —1, 0, 1 in the first 
Brillouin zone. It is noted that in Eq. (19) we have included the contribution of topological 
phase for FM particles with trigonal crystal symmetry (i.e., 2irS/3).. The low- lying energy 
level spectrum, which corresponds to the splittings of m-th excited state due to the resonant 
quantum coherence of the magnetization vector between energetically degenerate states, is 
found to depend on the parity of total spin of FM particle significantly. If S is an integer, the 
low-lying energy level spectrum is tie*™ — 2hAe t ™, and he 1 ™ + hAe*™, the latter being doubly 
degenerate. But if S is a half-integer, the low-lying energy level spectrum is he 1 ™ — hAe 1 ™, 
and he 1 ™ + 2hAe t ™, the former being doubly degenerate. This spin-parity effect is the result 
of phase interference between topologically distinct tunneling path. 

At the end of this section, we discuss the possible relevance to the experimental test for 
spin-parity effects in single-domain FM nanoparticles. First we discuss the thermodynamic 



behavior of this system at very low temperature T ~ T = TiAe^ 1 jk B . For FM particles with 
trigonal crystal symmetry at such a low temperature, the partition function of the ground 
state is found to be 

Z = exp (-ptieF) [exp (±2phAe^ 1 ) + 2 exp (^/3hAe^)] , (20) 

where upper sign corresponds to integer spins, lower sign corresponds to half-integer spins, 
and 4" = wi/2. Then the specific heat is c = — T {d 2 F/dT 2 ), with F = -k B T In Z. For the 
low temperature case, the result is 



tri\ 



c = 18k B (pHAer) 2 e W (±phAe o l 

V J [exp(±2^A4") + 2exp(^A4")] 2 

In Fig. 1, we plot the temperature dependence of the specific heat for integer and half- 
integer spins at very low temperature < T/T < 1. It is clearly shown that the specific 
heat for integer spins is much different from that for half-integer spins at sufficiently low 
temperatures. When the temperature is higher hAe$ % <C k B T < hui, the excited energy 
levels may give contribution to the partition function. Now the partition function is 



1 + (l - e-^ 1 ) (>iA4") 2 J (2g ie -^ l/2 ) 



(22) 



for both integer and half-integer spins, where Z = 3e /^i/ 2 / (l — e l 3 ^ 1 ^ is the partition 
function in the well calculated for k B T <C AU over the low-lying oscillatorlike states with 
e% = (m + 1/2) co u and = 3y/2(V/hS) ^/KjT 2 - h (x) = En=o (x/2) 2n / (n!) 2 is the mod- 
ified Bassel function, and q\ = (2 9 / 2 /3^J \ 1 / 2 S > 1. We define a characteristic temperature T 
that is solution of equation q 1 e~ huJl l kBT _ ^ temperature T = hui/2 lngx characterizes 
the crossover from thermally assisted tunneling to the ground-state tunneling. In Fig. 2, 
we plot the temperature dependence of the specific heat for integer and half-integer spins 
at high temperature 30 < T/T < 60. The result shows that the spin-parity effect will be 
lost at high temperatures. The specific heat for integer spins is almost the same as that for 
half-integer spins. 



9 



IV. MQC FOR TETRAGONAL AND HEXAGONAL CRYSTAL SYMMETRIES 



In this section, we will apply the method in Sec. Ill to study spin tunneling in FM 
particles with tetragonal and hexagonal crystal symmetry. For the tetragonal symmetry, 

E (0, 0) = K 1 cos 2 9 + K 2 sin 4 9 - K' 2 sin 4 9 cos (40) + E , (23) 

where K x 3> K 2 , K' 2 > 0. The energy minima of this system are at 9 = it/ 2, and = 
0, 7r/2, 7r, 37r/2, and other energy minima repeat the four states with period 2-7T. The 
problem can be mapped onto a problem of one-dimensional motion by integrating out the 
fluctuations of 9 about 7r/2, and for this case V ((f)) = 2 (K' 2 V/h) sin 2 (20). Now V (0) is 
periodic with period 7r/2, and there are four minima in the entire region 2ir. The periodic 
instanton configuration with an energy E > is sin 2 (20 p ) = 1 — k 2 sn 2 (u 2 t, k), where 
k = yj(nl - 1) /nj, u 2 = 2 5 / 2 (V/hS) yj K X K' 2 , and m = ^2K' 2 V/TiE > 1. The associated 
classical action is S p — W + 2E/3, with 

W = 2 1 / 2 S\ § \E (k) - (l - k 2 ) K (k)} . (24) 

The general formula Eq. (4) gives the low-lying energy shift of m-th excited state as 



hAe^ = (K X V) \H m +l)s( m - l i) exp (^X^S) , (25) 

with A = K' 2 /Ki. The periodic potential V (0) can be viewed as a superlattice with lattice 
constant ir/2 and total length 2ir, and the Bloch theorem then gives the energy level spectrum 
of m-th excited state e 1 ^ = (m + 1/2) uj 2 as E% = — 2As t £ l cos [(S + £) n/2], where £ = 
— 1,0,1,2 in the first Brillouin zone. It is easy to show that the low-lying energy level 
spectrum is he^ ± 2hAe t 1 ^, and he 1 ^ for integer spins, the latter being doubly degenerate. 
While the level spectrum is he^ ± V^hAe^ with doubly degenerate for half-integer spins. 
At a very low temperature T ~ T = hAe^/ks, the specific heat is 

c = 4k B (phAel e ) 2 - ]— - , . , (26a) 

v > 1 + cosh (2/3hAel e ) ' v ; 

for integer spins, while 
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for half-integer spins. The tunneling behavior for integer spins is almost the same as that 
for half- integer spins at high temperature HAe^ <C k B T < hu 2 . 
For the case of hexagonal symmetry, 

E (6, 0) = K x cos 2 9 + K 2 sin 4 9 + K 3 sin 6 9 - K' 3 sin 6 9 cos (60) + E , (27) 

where K\ 3> K 2 , K 3 , K 3 > 0. The easy directions are at 9 — n/2, and = 0, 7r/3, 27r/3, 7r, 
47r/3, 57r/3, and other energy minima repeat the six states with period 2ti. For the present 
case, V (0) = 2 (K' 3 V/K) sin 2 (30) is periodic with period 7r/3, and there are six minima 
in the entire region 2n. The periodic instanton configuration at a given energy E > is 
sin 2 (30 p ) = 1 - A; 2 sn 2 {u 3 t, k), where k = ^J{nj - 1) /nj, w 3 =(3x 2 3 / 2 ) (V/HS) s]K x K' 3 , 



and rii = y2K 3 V/hE > 1. Correspondingly, the classical action is = + 2E(3, with 

W = 2 3/2 1 SW ^ [E (A;) - (l - A; 2 ) K (jfc)l , (28) 
and the low-lying energy shift of m-th excited state is 



A:f ; ; = J-^^L^^AsKD^-^exp f-^A 1 / 2 ^ , (29) 
m! 3 m "27T2 V 3 / 



with A = K' 3 / K x . Now V (0) can be regarded as a one-dimensional superlattice with lattice 
constant n/3. By applying the Bloch theorem and the tight-binding approximation, we 
obtain the energy level spectrum of m-th excited state e 1 ^ = (m + 1/2) uo 3 as E^ = e 1 ^ — 
2Ae^ cos [(S + £) 7r/3], where £ = —2, —1, 0, 1, 2, 3. If 5 is an integer, the low-lying energy 
level spectrum is Ue 1 ^ ± 2hAe^, and Ue 1 ^ ± hAe^, the latter two levels being doubly 
degenerate. If S is a half-integer, the level spectrum is he 1 ^ ± y/SfiAe 1 ^, and he 1 ^, all three 
levels being doubly degenerate. Then the specific heat at sufficiently low temperatures is 

, ln2 [4 + 4 cosh (phAe^) + cosh ( 2(3hAe't e ) cosh ((3hAe^ e ) 

c = 2k B ((3hAe h *) ± V z ~r V ,° J -A (30a) 

[cosh (2f3hAe% e ) + 2 cosh [/3hAe^ 

for integer spins, while 
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, . X 2 2 + cosh (V3f3hAe% e ) 

c = 6k B (f3hAe h e ) ) ^T2, (30b) 

v ' [l + 2 cosh (v^ftAe^ 



for half-integer spins. 

In brief, the low-lying energy level spectrum and the heat capacity of the magnetic 
tunneling states for tetragonal and hexagonal symmetry are found to depend on the parity 
of the total spins, resulting from the Wess-Zumino phase interference between topologically 
distinct tunneling paths. And this spin-parity or topological phase interference effect will 
be lost at high temperatures. 



V. CONCLUSIONS 

In summary, we have investigated the topological phase interference effects in spin tunnel- 
ing at excited levels for single-domain FM particles with trigonal, tetragonal, and hexagonal 
crystal symmetries. The low-lying tunnel splittings between m-th degenerate excited states 
of neighboring wells are evaluated with the help of the periodic instanton method, and the 
energy level spectrum is obtained by applying the Bloch theorem and the tight-binding 
approximation in one-dimensional periodic potential. 

One important conclusion is that for all the three kinds of crystal symmetries, the low- 
lying energy level spectrum for integer total spins is significantly different from that for 
half-integer total spins, resulting from the Berry phase interference between topologically 
distinct tunneling paths. For FM particles with simple biaxial crystal symmetry, which 
has two degenerate easy directions in the basal plane (i.e., the double- well system), it has 
been theoretically shown that the tunnel splitting is suppressed to zero for half-integer spins 
due to the destructive phase interference between topologically different tunneling paths 
connecting the same initial and final states. However, the structure of low-lying tunneling 
level spectrum for the trigonal, tetragonal, or hexagonal crystal symmetry is found to be 
much more complex than that for the biaxial crystal symmetry. The low-lying energy level 
spectrum can be nonzero even if the total spin is a half-integer for the trigonal, tetragonal, or 
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hexagonal crystal symmetry. Note that these spin-parity effects are of topological origin, and 
therefore are independent of the magnitude of total spins of FM particles. The heat capacity 
of low-lying magnetic tunneling states is evaluated and is found to depend significantly on 
the parity of total spins for FM particles with different crystal symmetries at sufficiently low 
temperatures, providing a possible experimental method to examine the theoretical results 
on topological phase interference effects. And the spin-parity effects will be lost at high 
temperatures. Our results presented here should be useful for a quantitative understanding 
on the topological phase interference or spin-parity effects in resonant quantum tunneling 
of magnetization in single-domain FM particles with different crystal symmetries. 

More recently, Wernsdorfer and SessolS have measured the tunnel splittings in the 
molecular Fe§ clusters, and have found a clear oscillation of the tunnel splitting with the 
field along hard axis, which is a direct evidence of the role of the Berry phase in the spin 
dynamics of these molecules. It is noted that the theoretical results presented in this paper 
are based on the instanton method, which is semiclassical in nature, i.e., valid for large spins 
and in the continuum limit. Whether the instanton method can be applied in studying the 
spin dynamics in molecular clusters with S = 10 (such as Fe 8 ) is an open question. 

The theoretical calculations performed in this paper can be extended to the FM and 
AFM particles in a magnetic field. Work along this line is still in progress. We hope that 
the theoretical results presented here will stimulate more experiments whose aim is observing 
the topological phase interference effects in nanometer-scale single-domain magnets. 
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Figure Captions: 

Fig. 1 The temperature dependence of the specific heat for integer and half-integer spins 
at very low temperature < T/T < 1. 

Fig. 2 The temperature dependence of the specific heat for integer and half-integer spins 
at high temperature 30 < T/T Q < 60. Here A = 0.001. 
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